Quantum coherence in ion channels: Resonances, Transport and Verification 
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Recently it was demonstrated that long-lived quantum coherence exists during excitation energy 
transport in photosynthesis. It is a valid question up to which length, time and mass scales quantum 
coherence may extend, how to one may detect this coherence and what if any role it plays for the 
dynamics of the system. Here we suggest that the selectivity filter of ion channels may exhibit 
quantum coherence which might be relevant for the process of ion selectivity and conduction. We 
show that quantum resonances could provide an alternative approch to ultrafast 2D spectroscopy 
to probe these quantum coherences. We demonstrate that the emergence of resonances in the 
conduction of ion channels that are modulated periodicalUy by time dependent external electric 
fields can serve as signitures of quantum coherence in such a system. Assessments of experimental 
feasibility and specific paths towards the experimental realization of such experiments are presented. 
We show that this may be probed by direct 2-D spectroscopy or through the emergence of resonances 
in the conduction of ion channels that are modulated periodically by time dependent external electric 
fields. 



PACS numbers: 

A wide variety of transport processes on all length and 
time-scales are of fundamental importance for the func- 
tion of biological systems. For a long time it was as- 
sumed that these transport processes may be modelled 
accurately as classical random processes and the possible 
relevance of non-classical phenomena such as quantum 
coherence was mostly ignored. The experimental demon- 
stration of the presence of long-lived quantum coherence 
in the photosynthetic energy transfer at cryogenic tem- 
peratures [2 [5] and room temperature has forced a 
critical reevaluation of this position. As a consequence, 
the detailed experimental and theoretical study of quan- 
tum coherence in bio-molecular systems and its possible 
relevance to explain functional properties of those sys- 
tems is currently receiving rapidly increasing attention. 

The transport in the photosynthetic complex or the 
reactions during the first steps in vision [4l [5] involve ini- 
tiation and fast transfer of electronic excitation via an 
absorbed photon on the femtosecond time scale. Given 
that these systems have been evolutionary optimized to 
function at highest efficiencies (e.g. human vision has few 
photon sensitivity [SI [Z] ) , they might seem to naturally 
suggest themselves as areas to look for quantum coher- 
ence. More generally however, on time scales that are 
short compared to the decoherence time of an observed 
biological process quantum mechanical coherence can be 
expected to occur and play a role in biological function 
despite (or even because of) the strong interaction of bio- 
molecular systems with their environment. 

Indeed, there may be non-trivial effects, positive or 
negative, due to the interplay between quantum coherent 
dynamics and decoherence from environments. In models 
of photosynthetic energy transfer it has been shown that 
compared to the purely quantum coherent case the effi- 
ciency of excitation transfer is enhanced by the presence 



of dephasing noise |8l [9] originating from an environment 
of vibrational modes. It was shown that the transport 
performance of the system is optimal when it is neither 
fully decohered and hence classical nor fully quantum 
coherent. The mechanisms underlying this phenomenon 
have been identified and understood to apply more gen- 
erally beyond the specific setting of energy transport in 
photosynthetic complexes [lOL and optimized trans- 
port settings in simpler structures have been determined 
|12) . The role of quantum coherence in other light har- 
vesting systems [T3] and for speed-up at early stages of 
the transport process has been examined [14j and the 
presence of entanglement [lOl [15] has been studied in 
some detail. 

These results provide further motivation to explore 
the relevance of quantum coherence and decoherence for 
transport in other biological entities. Here we discuss 
another system in which the interplay between quantum 
coherence and environmental noise could be responsible 
to explain some of the functional observations of the sys- 
tem. 

Ion-channels are protein complexes that regulate the 
flow of particular ions across the cell membrane and are 
essential for a large range of cellular functions [TB]. Be- 
sides their role in neuronal communications in which volt- 
age gated channels and ligand gated channels are involved 
in the generation of action potentials and mediating 
synaptic release, more generally ion-channels play a key 
role in processes that rely on fast responses on the bio- 
molecular scale. Examples include muscle contraction, 
epithelial transport and T-cell activation [THl [TT] . Struc- 
turally, ion-channels are integral transmembrane protein 
complexes with multiple subunits whose relative spatial 
arrangement forms a pore through which ions can flow 
in or out of the cell along a concentration gradient. A 
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common feature in a large number of ion-channels is the 
presence of a gate that can be activated by chemicals, 
voltage, light or mechanical stress and a selectivity filter 
which is responsible for allowing only a specific type of 
ions to pass. The selectivity filter is only a few angstroms 
wide hence the ions have to move through it in a single 
file fashion without their hydration shell jl8J . 

A large number of ion-channels have been success- 
fully crystallized over the last decade and x-ray crys- 
tallographic data have provided us with structural un- 
derstanding of these protein complexes on the atomic 
level . In parallel based on this data and by using 

molecular dynamics simulation different functional mod- 
els for the mechanisms of ion selectivity and conduction 
have been developed [25U27j . however experimental ob- 
servation of these dynamics in real time and on atomic 
level has remained challenging. 

Amongst the best studied ion-channels are potassium 
channels. Their family includes a wide range of channels 
with different conductivities and gating mechanism in- 
cluding the voltage gated potassium channels which are 
responsible for restoring the membrane potential during 
the course of action potentials in neurons |16j . While 
the gating mechanism can greatly vary across different 
types of potassium channels, the sequence of amino acids 
forming the selectivity filter is fairly conserved across all 
potassium channels. The structure of the selectivity fil- 
ter has been best studied in the bacterial KcsA channel 
[TO] on which the following discussions are focused. 

During its open gate state an efflux of 10® ions per sec- 
ond can be experimentally observed |28j which is close to 
the diffusion limited rate. At the same time the selec- 
tivity filter is capable of selecting potassium over sodium 
with a ratio of at least 10"* ; 1 [19 . Furthermore, it is 
known from the x-ray crystallographic data of the KcsA 
that its selectivity filter forms a pore with the diameter 
of ^ 0.3 nm and a length of ^ 1.2 nm. This implies that 
individual potassium ions cannot move through the se- 
lectivity filter with their 8 water molecule hydrated shell 
which must be shed before entering the selectivity filter. 
The transfer through the filter then proceeds in single 
file fashion where potassium and water molecules alter- 
nate. This makes the mentioned observations on trans- 
mission and discrimination rates even more remarkable. 
In a number of models of the selectivity filter this fact 
is accounted for by assuming that the ion transport is 
achieved via an interplay of ion - filter attraction and ion 
- ion repulsion while the selectivity is achieved through a 
mechanism which exploits the lower water shell dehydra- 
tion energy for potassium than for sodium [23, 28, 29J. 
Certainly, the magnitude of the thermal fiuctuations of 
the backbone atoms forming the selectivity filter is large 
relative to the small size difference between Na and K, 
raising fundamental questions about the mechanism that 
gives rise to ion selectivity. This suggests that the tradi- 
tional explanation of ionic selectivity should be reexam- 
ined [H]. 

A closer look at the involved dimensions and energet- 




FIG. 1: Schematic illustration of the KcsA postassium chan- 
nel after PDB 1K4C. KcsA protein complex with four trans- 
membrane subunits (left) and the selectivity with four axial 
trapping sites formed by the carbonyl oxygen atoms in which 
a potassium ion or a water molecule can be trapped. 



ics of the process reveals that the underlying mechanism 
for ion transmission and selectivity might be not entirely 
classical. Fig{l] shows the structure of the KcsA and the 
details of the selectivity filter. It is formed by four chains 
of amino acids each from one of the four protein com- 
plex subunits constituting the ion channel. Each chain 
is made of five amino acids from whom a carbonyl group 
oxygen atom is pointing toward the pore (Figjl] inset). 
This configuration leads to a series of axial potential min- 
ima (i.e. binding sites) and maxima for the potassium 
ions. In each of these binding sites which are formed by 
four in plane oxygen atoms either a dehydrated potas- 
sium ion or a water molecule can be trapped (Figjljinset). 
Numerical calculations of free energy maps of the ion con- 
duction shows two energetically nearly degenerate config- 
uration paths for the ion transport which are separated 
by only ^ 2kcal/mol [23]. The transport along each con- 
figuration path involves a coordinated dislocation of wa- 
ter -potassium chain. During these transitions different 
sequences of water and potassium occupy the four bind- 
ing sites. This is also supported by the x-ray diffraction 
data which has revealed two equally represented popula- 
tions corresponding to two different sequences of water 
and potassium in the filter which are referred to as 1,3 
and 2,4 state [55]. In this arrangement the axial separa- 
tion of the potential minima is ^ 0.24 nm and the height 
of the barrier fluctuates based on the presence or absence 
of an ion at a particular site and the thermal vibrations of 
the protein in the range of ~ 1 — 5 kcal/mol correspond- 
ing to ~ 1.7 — SksT [27j. This situation allows for at 
least two different but equivalent considerations through 
which quantum coherence could be in involved in the ion 
transmission process. 

First, based on the thermal energy oi E — kBT/2 w 
2 X 10^^^ J of the ions the de-Broglie wavelength Xdb = 
h/\/2mE of the potassium matter wave can be estimated 
to be ^0.05 nm. Given that the periodicity of the axial 
modulation of the potential is '-^0.25nm it is within a 
factor of five of the potassium matter wavelength. In 
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this sense the transmission could involve a process which 
can be viewed as the diffraction of the potassium matter 
wave off a one dimensional axial "grating" formed by the 
modulations of the potential energy. 

Alternatively the transmission could be seen to in- 
volve quantum tunneling through the potential barrier 
between individual neighboring binding sites. This can 
seen by estimating the tunneling probability, ptun = 
^-A^27n{Eo-E)/h^ ^ ^^[f^l^ ^]^g energy of a potassium ion, 
TO, the mass of potassium, Eq — E the height and A the 
width of the tunneling barrier. Using the parameters of 
the selectivity filter [211 HZ] this yields a tunneling prob- 
ability of ~ 7 X 10~^ for a single "try". Multiplying this 
value with the trapping frequency of a bound ion in a 
harmonic potential of ~ lO^^Hz yields an overall trans- 
mission approaching the observed rates of 10^ ions per 
second. This estimate assumes a square well potential. 
However, a more accurate description would assume a 
quartic double well potential V{x) = ax'^ — 2Px^ whose 
parameters are chosen such that the potential barrier at 
V{0) has a height AE = OMeV = l.bkT, and a sep- 
aration of the minima of ^ 0.24nTO. Furthermore, as- 
suming that the degree of excitation inside a potential 
well will be compatible with the thermal energy, we can 
then compute the energy level splitting between levels 
that are split due to tunneling. This splitting gives the 
tunneling rate for a system in that energy level. Typical 
values range between '-^lO'^s^^ and ^10"'^°s~"'^ for energy 
levels around the thermal energy. In these estimations 
we have assumed a static potential and have ignored any 
rotational vibration modes in the carbonyl groups. How- 
ever, it has been theoretically shown [571 130] that these 
could significantly lower the potential barriers, so that 
the tunneling probabilities given above would this sense 
represent only a lower bound. 

Based on these estimations we hypothesize that quan- 
tum coherence could be present in the selectivity filter 
and might be involved in explaining some of the func- 
tional features such as the high ion conduction rate and 
ion discrimination rate. However, given the strong inter- 
action of the system with the environment these coher- 
ences may be short lived. A rough estimate suggests that 
the decoherence time will be on the order of 10~^ — 10~®s 
and thus of the order of the expected tunneling rates [31] . 
In fact, we expect that the interplay between quantum 
coherent dynamics and decoherence from the environ- 
ment might actually be necessary for explaining the dy- 
namics of the selectivity filter. As demonstrated for the 
excitation transport through photosynthetic complexes 
the presence of dephasing noise [8, 9 may even enhance 
the efficiency of transport compared to the purely quan- 
tum coherent case. 

A key feature of such coupled quantum systems that we 
will be exploiting here is that interference, constructive 
or destructive, in the underlying dynamics may lead to 
resonances in the transport cfhciency. These resonances 
allow us to infer the presence of quantum coherence as 
they are absent in pure rate equation models. These res- 



onances may be measured directly in the energy, particle 
or charge current that the systems allows - the measure- 
ment of correlations in time or coherence in relation to 
the driving field are not required. 

This represents two distinct advantage. First the direct 
detection and verification of quantum coherence in bio- 
logical systems has been so far been based mainly on two 
(or higher) dimensional femtosecond spectroscopy tech- 
niques 1,32, in which quantum coherence manifests itself 
as a beating signal in the intensity of the cross peaks. In 
two-dimensional femtosecond spectroscopy three pulses 
and a strongly attenuated local oscillator are incident 
on the sample. The first pulse creates a coherence that 
evolves for a time period usually referred to as the evolu- 
tion or coherence time, then the second pulse creates an 
excited-state population that evolves over a second time 
period called the waiting time or mixing period. Finally, 
the third pulse generates a coherence that accumulates 
phase in the opposite direction. Thereafter re-phasing 
occurs and a signal pulse is emitted in a direction deter- 
mined by the phase matching condition. The coherence 
time and the detection time are scanned in a range of 
a few hundred picoseconds in the vibrational 2D spec- 
troscopy and in the range of femtoseconds for studies in 
the visible range while the signal is measured through 
a heterodyne detection for each combination of the co- 
herence and the detection time. In this way coupling 
strengths and energy transfer in molecules are probed. In 
the absence of coupling contributions from the excited- 
state absorption and emission cancel each other and as a 
result no off-diagonal peaks appear in the spectrum. In 
the presence of coupling the contributions from the ex- 
cited state absorption and emission do not fully cancel 
and a so-called "cross peak" emerges in the spectrum. 
Intuitively speaking, the 2D spectroscopy can probe the 
"memory" of a system. As a result it can be used to ex- 
amine how fast an initially coherent state decoheres with 
time. 

Even though several laboratories are now mastering 
this method in the visible [Iti3l [33] and in the IR regime 
[52H57| its implementation remains an experimental chal- 
lenge. Besides the usual complexity involved in the de- 
sign, operation and maintenance of a coherent ultra- 
fast 2D spectroscopy system, it implies stringent require- 
ments for stability and precision for the inter-pulse delay 
times. 

Second, in some systems, such as in current-carrying 
polymers or in the selectivity filter the observation of 
resonances in transport rates may not only provide an 
alternative approach towards the verification of the exis- 
tence of quantum coherence but also a means for direct 
demonstration of its importance for the system dynamics. 
Indeed, a further potentially significant implication of our 
work is that the presence of quantum coherence for exam- 
ple in the selectivity filter and the concomitant existence 
of sharp resonances in the transport may point towards 
a novel mechanism for the explanation of the high degree 
of selectivity. Indeed, relatively small changes in the sys- 
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tern parameters (mass, ion radius, interaction strengths 
of ions traversing the channel) may then lead to sharp 
changes in the conductance. 

Quantum coherence & transport - The dynamics of ion 
conduction inside an ion channel is highly complex due to 
the manifold interactions between the ions and the mul- 
tiple degrees of freedom of the channel. In the following 
quantum mechanical analysis of a transport process it is 
our aim to avoid this complexity while retaining some 
structural elements of the dynamics of an ion channel. 
This will allow us to show that in this simplified model 
quantum coherence can give rise to very significant effects 
in the conductance of the channel. In turn this suggests 
that a more realistic model may also exhibit similar fea- 
tures. This provides the motivation for considering pos- 
sible experimental set-ups that would be able to verify 
or falsify the hypothesis of the presence of coherence in 
a real ion channel. 

The theoretical and numerical discussion, first defines 
the basic dynamical equations governing the transport 
and discusses their relationship to real systems. Then we 
move on to demonstrate that coherently driven transport 
channels exhibit resonances in their conductance which, 
in turn, can be linked to the level of quantum coherence 
in the system. In fact, we will elucidate the mechanism 
behind these resonances and hence understand why they 
do not occur in systems exhibiting classical rate equation 
dynamics. Based on the insights that we gain from this 
analysis we then proceed to discuss in some detail the ex- 
perimental feasibility of our predictions and discuss the 
concrete settings underlying our current experimental ef- 
forts in this direction. 

The basic dynamical equations ~ In the present con- 
text we are focusing on quantum dynamical equations 
that are motivated by the example of the selectivity fil- 
ter of ion channels. Nevertheless, their essential features 
should capture a wide variety of transport processes in 
bio-molecular systems. One of the underlying assump- 
tion is that to a good approximation the transport is 
particle number conserving, that is, it is described by 
a tight-binding type Hamiltonian that is mainly subject 
to dephasing-type noise. This is a reasonable model as- 
sumption for excitation transport in photosynthesis (ex- 
citation annihilation rates are short compared to the ob- 
served transport times) or the transport of ions through 
an ion channel (as ions rarely escape sideways through 
the channel walls). As we are principally motivated by 
ion channels we restrict attention to a linear chain of sites 
in which an excitation may be exchanged between near- 
est neighbors even though the ideas presented here are 
not limited to that setting. Note that an excitation at a 
particular site in the following Hamiltonian model does 
not need to represent the presence of a potassium ion at 
a certain position in space. Rather each site may also 
be thought of as specific potassium-water single file con- 
figurations that may arise during the transport. Each of 
these configurations will sit in an effective potential well 
and is separated from other configurations by a potential 



barrier that it may traverse through quantum coherent 
tunneling or thermal activation |18j . The coherent part 
of the evolution is then described by the Hamiltonian 
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with hopping rates Cfc(t) and site energies Huikit) that 
may vary across sites and in time. Let us furthermore 
assume that the first and last site of the chain are con- 
nected to environments that insert or remove excitations 
from the system or, in other words, drive transitions be- 
tween configurations. The first site may for example be 
linked to a source that supplies particles (electrons from 
a lead or potassium ions from the cytoplasm entering the 
ion channel thus driving the system to a different con- 
figuration) while the last site may then be linked to a 
similar environment but with lower voltage and/or lower 
concentration of particles. In general one may describe 
such processes by Lindblad [SI] operators of the type 

+^s/dns/d[-<^I<^tp - Pt^if^t + ScTj^pcrf] (2) 

where s (d) stand for source (drain), (F^j) for the cou- 
pling strength between the chain and the environments 
and where Us {rid) represent concentrations. Terms of the 
form ai p — pa^a'^ are responsible for the evolution 
of the system when no excitation has entered or left the 
system, while the term pa^ {2a^ pa^) describes the 
event when an excitation has left (entered) the system. 

Dephasing noise also originates from an interaction 
with an environment, for example with the oscillatory 
degree of freedom of the carbonyl groups in the ions 
channel, but does not create or destroy excitations in the 
chain. Dephasing can be described by another Lindblad 
operator of the form 
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with site dependent dephasing rates jk- It should be 
noted that this model ignores temporal and spatial cor- 
relations in the noise. While such correlated noise may 
lead to quantitative changes in the following considera- 
tions the basic principle of using the presence or absence 
of resonances in the conductivity to determine the pres- 
ence of quantum coherence remains valid. Besides the co- 
herent hopping between the sites there may also be ther- 
mally activated transitions between nearest neighbours. 
These would again be described by Lindblad operators 
of the form 

N~l 
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FIG. 2: A source drives excitations into the first site of a chain 
where nearest neighbors are coupled by a hopping interaction. 
The excitations leave the chain via the last site. 



Hence the most general form of quantum dynamical equa- 
tion that we will be considering here is of the form 

^ = -'-[H,p]+Csip)+Cdip)+Cth{p)+Cdeph{p). (5) 

To determine the conductivity of such a system we will 
consider Psink{t) which may be obtained directly by in- 
tegrating the population of the last site multiplied with 
the net transfer from that site into the sink, that is 
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'^^dPn,n{t)dt. 



As the current we will then define the asymptotic rate of 
growth of Psink (t) that is 



lim ^:P^(t). 



(6) 



Quantum coherence & conductivity - Now we would 
like to demonstrate that the conductivity of an externally 
driven channel exhibits resonances at which the conduc- 
tivity is strongly suppressed. We will then show that the 
depth of these resonance reduces in the presence of de- 
phasing noise and that it can be related to the amount 
of quantum coherence in the channel. Thus the depth of 
resonances may be taken as a measure for the presence 
of quantum coherence. 

We consider a channel subject to a constant and a 
time dependent potential due to applied electric fields. 
Furthermore we assume, for simplicity and clarity of the 
argument, that the inter-site hopping rates are constant 
across the chain, i.e. Ck = c. In this case the coherent 
part of the dynamics is given by a Hamiltonian of the 
form 
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which in turn implies that c^ /fto ^ 10^s~^ if c ^ f^o as 
is the case for an ion channel. 

In the limit of very long, ideally infinite chains, this 
Hamiltonian is known to exhibit the quantum coherent 
phenomenon of dynamic localization or dynamical sup- 
pression of tunneling, [35j |36]. The existence of this 
effect and its coherent character may be inferred from 
the following arguments. Consider a transformed pic- 
ture in which the time-dependent on-site energies van- 
ish at the expense of time dependent coupling strengths 
between the neighboring sites [55]. This is achieved 
by defining A{t) = —ftot — {^i / ui) sm{uit) to move to 
an interaction picture via the transformation \4^(t)) = 
g-iA(t) J2k ^'^t'^k \tl){t)). In this interaction picture the dy- 
namics is governed by a Hamiltonian 
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Here MIq is the energy difference between adjacent sites. 
In the ion channel this would be related to the energy 
difference between two close lying energy levels in adja- 
cent potential wells. This difference may not exceed the 
level spacing in one of the potential wells which in turn 
tends to be below 10^^ s~^. To achieve a transfer rate 
of the order of lO^s""'^ through the channel, the effective 
hopping rate between sites should be of the same order 
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where now the coupling rates are time-dependent. For 
small c the hopping dynamics between neighbouring sites 
is slow compared to the time dependence e**'^'^*^ and av- 
eraging this Hamiltonian over the interval [— tt/w, tt/o;] 
we find that the effective Hamiltonian takes the form 

JV-i ni _ _ 

Hi/h= cJna{-^){(T+a-^^^+a-f^a+^-^). (9) 

— ^ uJ 

/c=l 

Hence we expect that if flo — nuj for some n e M and 
Ei/uj coincides with a zero of the Bessel function J„, 
then the evolution of a wave-packet becomes periodic in 
time and the spreading of the wave-packet and hence 
transport is suppressed. For an infinite chain this can be 
demonstrated rigorously [351 ISS] ■ 

For a finite system that is in contact with a source 
and a sink, excitation transport will not vanish ex- 
actly. Firstly, within a time interval [— vr/o;, tt/o;] the 
wavepacket will not be stationary but will oscillate chang- 
ing its width. The localization length (the maximal ex- 
tent over which the wavepacket may spread in that in- 
terval) may exceed the length of the channel and hence 
leave it into the sink. On the other hand the quantum 
coherence responsible for the localization will be per- 
turbed by the incoherent processes that are taking place 
at source and sink. Nevertheless, it can still be expected 
that transport is suppressed significantly, even for short 
chains. Indeed, a numerical analysis confirms the pres- 
ence of dynamic localization even for the dynamics of a 
short chain. Fig. |3] indeed demonstrates that with just 5 
sites (similar results may be observed for shorter chains 
accompanied by higher oscillation frequencies w), high 
contrast resonances in the conductivity can be observed 
and that these coincide very well with the zeros of the 
relevant Bessel-functions. 

Note that the effective coupling rates in eq. ^ 
emerge through averaging over transition amplitudes and 
it hence becomes clear that the suppression of trans- 
port is a coherence effect. Another way to understand 
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FIG. 3: Conductivity versus strength Ei of the apphed ac- 
potential for a chain with parameters N = 5, Qo = 256 • 
W^s-\ = r, = W^s-\ LU = 2 - W^s-\ TIT = 1 and 
c = 8 ■ 10*s~^. For increasing dephasing rates 7 — 0, 5 • 
10^s-\l.5-10**s-^ and2-10*s"^ one observes a reduction in 
the depth of conduction resonances. 



these results is from the observation that the ac-field 
will let the onsite energies oscillate so that periodically 
in time the neighboring sites become energetically close 
enough so that a coupling c that is insufficient to bridge 
an energy-gap of the order of Eq will then induce tran- 
sitions. In fact, this is the periodic analogue of the 
mechanism by which dephasing helps to overcome en- 
ergy gaps [TU]. However, in the periodic case, the wrong 
timing, can lead to coherent suppression of transport as 
well 0ni- The physical origin of this trapping, and its 
quantum nature, may be understood by considering the 
extreme case of the Hamiltonian H/h = J2k=i f [1 + 
sign{sin{ujt + kTT))]{a'^ a'jl_^^ + tr^Tcr^^.]^) that is the cou- 
plings are switched on and off with a frequency w alter- 
natingly between neighboring pairs of sites. Within this 
model it is now easy to see why time-varying coupling 
rates may suppress transport in a quantum setting. For 
a; = c we observe that exp~'^=^^ '^/^ |1) = — 11) and as a 
consequence excitations do not propagate at all. Hence, 
we expect that for certain frequencies transport will be 
suppressed 

All these arguments show that the dynamic suppres- 
sion of transport relies on the existence quantum coher- 
ence, either through the destructive interference of tran- 
sition amplitudes or the existence of full oscillations be- 
tween sites. In purely classical rate equation models or 
strongly decohered quantum systems these effects will 
not exist. There is no destructive interference nor can 
the population between neighboring sites oscillate per- 
fectly as it will rather approach an equilibrium value in 
which both sites will be populated. As a consequence we 
will not observe complete suppression of transport in a 



classical rate equation models with time-modulated cou- 
pling rates between sites. In this sense dephasing may 
in fact assist the transport as it suppressed destructive 
interference at the resonance points. 

The expected sensitivity of this phenomenon to the 
presence of decoherence can be seen in Fig. |3] Decreasing 
depth of the conductance resonances is associated with 
increasing levels of decoherence which in turn can be ex- 
pected to be associated with decreasing levels of quantum 
coherence in the system. This is indeed a crucial obser- 
vation as it allows us to turn the situation on its head 
and use the depth of the contrast of the resonances as a 
tool to determine quantitatively the presence or absence 
of quantum coherence in the system from measurement 
of its conductivity. 
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FIG. 4: Incoherence versus conductivity at the first resonance 
of fig. [3] The best fit to the data suggest that incoherence is 
proportional to conductivity. 

There is no unique measure for the level of quan- 
tum coherence or incoherence in a state but the ob- 
served direct relationship hold true for a variety of mea- 
sures. As a measure of incoherence, the quantity C ~ 
J2k^i \Pk.kPi,i- Pk,iPi,k\ has the desirable properties that 
it vanishes for pure states, increases with decreasing size 
of the diagonal elements and is maximized for maximally 
mixed state. In addition C is concave, that is, it increases 
under mixing of quantum states. For the first resonance 
in Fig. [Sj comparing the depth of the conductivity res- 
onance with the depth of the resonance in the quantum 
incoherence C for different levels of dephasing noise, we 
find Fig. [4j This result is very well fitted by C ~ /. It 
is this relationship that we would like to exploit for the 
detection of quantum coherence from the observation of 
currents in a modulated ion channel. 

Theory Discussion & Speculation - The results above 
demonstrate that systems that are governed by coherent 
quantum dynamics will, for certain choices of parameters. 
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exhibit resonances in their transport efRciency, i.e. con- 
ductivity. As these resonances are absent in purely clas- 
sical rate equation models or in quantum systems subject 
to very strong decoherence, the observation of currents, 
for example in periodically driven ion channels, provides 
an interesting alternative for the observation and verifi- 
cation of the presence of quantum coherence. 

These resonances, caused by quantum coherence, are 
relatively sharp features in parameter space. As a con- 
sequence, in the presence of quantum coherence, the dy- 
namics is much more sensitive to small variations of sys- 
tem parameters such as ion mass, ion radius etc. In 
purely classical models the rather small differences in di- 
ameter between potassium and sodium ions are not ex- 
pected to result in significant changes in the ion conduc- 
tion. In addition to the effects on the conduction rate pre- 
sented above one may speculate, whether the enhanced 
sensitivity of the quantum dynamics may contribute to 
the selectivity of the filter of the ion channel. 

Needless to say, the dynamics of ion conduction in ion 
channels is in fact a many-body phenomenon in which 
many degrees of freedom interact. Recent numerical 
simulations of the ion-transfer j^T] suggest for example 
that the multi water-potassium chain does not move as 
one unit but rather in discrete units of individual water- 
potassium translocations, in which the translocation of a 
potassium ion to the next site is followed by the water 
molecule in picosecond time scale. Given the short time 
scale for this substructure of the transport process it pro- 
vides an additional mechanism through which quantum 
coherence may appear. 

Of course such dynamics is considerably more com- 
plex than the model Hamiltonians described here more 
rigorous conclusions in that direction need to await more 
detailed simulations and experimental investigations |30| . 

Experimental realization - As mentioned earlier, it 
is an advantage of the model presented above that it 
provides experimentally testable predictions for the ex- 
istence of coherent excitations in the selectivity filter 
through the observation of driven quantum resonances. 
This is experimentally particularly attractive as the pres- 
ence of coherence can be inferred via the measurement 
of bulk signals, eliminating some of the challenges asso- 
ciated with alternative and more direct methods such as 
2D femtosecond spectroscopy. As shown in FigjS] in the 
presence of coherence at certain frequencies of an external 
driving field, the ion conduction through the selectivity 
filter can be highly impaired. This would be in sharp 
contrast to what should be observed when dephasing is 
destroying all coherence and the conduction of ions is 
governed by pure rate processes making it essentially in- 
dependent of the amplitude of the alternating external 
driving force. 

There might be several ways for testing this predic- 
tion experimentally. In the following we would like to 
discuss one of these in some detail to demonstrate the 
in-principle feasibility of such measurements. One could 
for example be using state of the art electrophysiology 



methods with some modifications for adapting them to 
the particular conditions that are required for resonances 
in ion channels to be observed. Electrophysiology tech- 
niques such as patch clamping 42 have found applica- 
tions in a large number of biophysical studies for record- 
ing single ion channel currents. During the open gating 
period of a few ms these currents can be observed ex- 
perimentally in different patch clamp techniques such as 
the cell attached mode and the excised patch configu- 
ration [42]. In these methods a glass pipette with an 
inner diameter of lfj,m is filled with an electrolyte and 
is used to make a high resistance (> lO^fi) seal with 
the cell membrane. Generally, only a few ion channels 
can be present on a membrane patch of that size. Be- 
cause of the high resistance of the seal any exchange of 
ions between the electrolyte in the pipette and the bath 
(or the cell cytoplasm) passes through the ion channels. 
Voltage-gated channels have single-channel conductances 
typically range from 10 — 100p5. They last over the few 
ms open period of the gate are amplified by a current to 
voltage converter and recorded. 

Here we should mention that the fact that multiple 
ion channels could be present on a membrane patch of 
the size of pipette tip by no means affects the exper- 
imental signature of quantum coherence in a negative 
fashion. As discussed above the manifestation of quatum 
coherence in the selectivity filter would be observable as 
a change in the "bulk" ionic current through the channel 
as a function of the externally applied AC driving force. 
In that sense as long as the driving force is constant over 
the number of channels, it actually represents an advan- 
tage to have more channels within the patched area, as 
it would help in increasing the signal to noise ratio. 

But in order to be able to apply electrophysiology 
methods to experimentally drive the resonances in the 
selectivity filter and record the transmitted ion rate sev- 
eral conditions have to be met. These will be discussed 
in some detail in the remainder of this work. 

In contrast to many single ion channel studies in which 
investigators are interested in the gating dynamics of the 
channel, our recordings have to be made during the open 
state of the channel. This requires the reliable prepara- 
tion of the channel in its open gate state. Since the struc- 
ture of the selectivity filter is highly conserved across the 
potassium channel family the experiments can also be 
done in Kv voltage gated potassium channels 43j. For 
this category of potassium channels a number of pharma- 
cological agents are available which can be used to initi- 
ate the open state of the channel. Alternatively, a class 
of potassium channels, the calcium activated potassium 
channels [^5Hl5] can be used in which the open confor- 
mation can be trigged by adding calcium to the bath. 
Finally, it is also possible to introduce site mutations in 
the voltage sensor region or use enzymes which impair the 
function of the voltage sensor and lead to permanently 
open channels [46] . 

Next, it is necessary to establish an ionic current across 
the channel in the absence of the AC driving force. 
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Again, while this could be done through a current in- 
jection via the patch clamp amplifier, it would be more 
convenient and also reassemble the natural conditions by 
establishing a potassium concentration gradient across 
the ion channel. The role of the potassium channel dur- 
ing the action potential is to restore the membrane po- 
tential from ^ -t-SOrnV^ after the fast opening of the 
sodium channels back to the resting membrane poten- 
tial of ~ —70mV. During this period potassium flows 
out of the cell along its concentration gradient ( 90mM 
intracellular vs. 3mM extracellular). Given the equi- 
librium potential 16J of potassium of — 85mV at this 
concentration difference this represent a driving force of 
llbmV . For our experimental considerations it would 
not be necessary to hold the membrane potential at 
+30™!^. Once the potassium concentration difference 
( 90mM intracellular vs. 3mM extracellular) is gen- 
erated by adding more potassium to the patch pipette 
than to the bath, even at OmV membrane potential, a 
driving force of 85mV for the ions would be available. If 
we assume an excised patch configuration in which a the 
current through a few channels on a dissociate patch of 
membrane sealed to the pipette is measured, this driving 
force will establish a constant current of a few pA. This 
constant driving force and the associated ionic current 
represent Eq in equation ([T]). 

In order to observe the predicted quantum resonances 
via the ionic currents through the channel an alternating 
driving potential, Ei cos ujt with the frequency uj has to 
be applied across the channel. This frequency w should 
be chosen to exceed the effective hopping rate in our the- 
oretical model which in turn is expected to be of the order 
of the ion channel transmission rate, i.e. around lO^s^^. 
As argued before, for an energy difference IiQq between 
adjacent potential wells, the resonances will be found at 
the zeros of the Bessel function Jjio/tj- We can expect 
Hq ^ Lu and hence the first zero of the Bessel function 
Jno/uj is, in leading order, located at f^i/w = fio/w [47] . 
The required ac-field could be applied across the chan- 
nel by using an external function generator synchronized 
with the amplifier. Alternatively an RF field at the same 
frequency with its polarization aligned along the axis of 
the ionic propagation might be used as an alternating 
driving force. 

However, because of the high frequency of the AC field 
the membrane capacitance and the stray capacitance of 
the pipette glass will lead to additional capacitive cur- 
rents significantly higher that the membrane current. 
Moreover, the time constant r for establishing a potential 
difference across the membrane is given by r = CmRp, 
with Cm the membrane capacitance and Rp the pipette 
resistance. Therefore it would be desirable to minimize 
the area of the membrane patch across which the AC 
field is applied. This limitation excludes any whole cell 
recordings and limits our considerations to the excited 
patch configuration for which the membrane capacitance 
is minimized as it will be determined by the surface area 
of the pipette tip. (see Fig [5]). The excised patch config- 




FIG. 5: Equivalent cuircuit for the excised patch configura- 
tion and external driving force to induce quantum resonances. 

uration is also advantageous as it allows the study of the 
ion channels in isolation and provides better control for 
the application agents on both sides of the channel. 

Assuming an inner diameter of l^m and a typical 
membrane capacitance of Cm ~ l^F/crn^ this results 
in a membrane capacitance of Cm ^ 8 x 10^ ^^F. As we 
would like to apply an AC field at lOOMHz the time 
constant r should be at least 5ns. Based on the mem- 
brane capacitance of Cm ~ 8 x 10~^^F estimated above, 
it means that the pipette resistance Rp should not be 
larger than O.GMfl. Typical pipette resistances are in 
the range of tens of Mil, but it should be possible to re- 
duce these values by appropriate shaping of the pipette 
tip. 

Further, the mentioned capacitive currents required 
some careful considerations. As shown in Fig[5]the mem- 
brane capacitance Cm and the pipette capacitance Cp are 
both contributing to the induced capacitive currents by 
the fast alternating AC field. While the membrane capac- 
itance can be minimized by using the excised patch con- 
figuration, minimizing the pipette capacitance requires 
using special type of glasses [48], minimizing the level 
of the electrolyte in the pipette and coating the pipette 
with insulating raisins which would thicken the external 
pipette walls and reduce the stray pipette capacitance. 
This can be done by coating the pipette with a layer 
of Sylgard, an inert, hydrophobic, translucent elastomer 
raisin which is cured rapidly by heat [49) . These measures 
could reduce the pipette capacitance to Cp ~ 10~^^F. 
However, this would be still about an order of magni- 
tude larger than Cm and therefore the major contributor 
to the capacitive current. 

As discussed above the applied AC field has to induce 
a membrane potential change of 315mV in lOns. With 
Cp as the dominating factor this means that the AC field 
would lead to a capacitive current of 3fiA which is six 
orders of magnitude higher that the pA ionic currents 
through the channel. It could be somewhat challenging 
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to separate these two components, however the difference 
in the timing of these currents could provide the means to 
separate them. The above choice for the time constant r 
is such that the peak the peaks of the capacitive currents 
are separated by the period of the AC field. Therefore 
a lock in detection technique can be used to reject the 
capacitive currents and isolate the ionic current. More- 
over for a fixed frequency (e.g. lOOMHz) it is possible 
to design additional electronic circuits which introduce 
negative capacitance to compensate the capacitive cur- 
rents. Finally a number spectral analysis technique can 
be applied to isolate the ionic current from the capacitive 
currents. Another approach which could eliminate the 
pipette capacitance altogether is to use chambers with 
suspended planar phospholipids bilayers in which the ion 
channels have been reconstituted ^50, . In this technique 
the two macroscopic chambers filled with electrolytes are 
separated by an artificial lipid layer in which the ion chan- 
nels have been reconstituted. Ions can move along their 
concentration gradient only through the ion channels and 
are recoded and amplified with a same amplifier system 
as discussed above. The AC field can be applied between 
the two chambers, while the current through the mem- 
brane is monitored. Since in this approach the pipette is 



eliminated the only source for capacitive currents is the 
lipid membrane. However in such a setup the membrane 
area is typically larger than the in the patch pipette so 
that it might lead only to a moderate reduction of the 
capacitive currents. 

Finally, the recent development of NV-centers point 
to an alternative method for recording pico- Ampere cur- 
rents in biological systems with high spatial and temporal 
control [5Tj . 

On the positive side the fact that the flux rate of ions 
at a given concentration is fairly constant and the fact 
that the width of the expected quantum resonances in 
Fig. |3]is not very sharp, is expected to help finding these 
resonance experimentally. 

Overall the above assessments show that currently 
available techniques from single channel biophysics stud- 
ies with moderate modification can be adopted and used 
for experimental investigations of quantum resonances in 
the selectivity filter of potassium channels. 
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To estimate the expected decoherence time we assume 
that the dominant form of noise stems from the stretch 
mode of the carbonyl groups in the ion channel. Their res- 
onance frequency is approximately ujco ~ 3.2 • 10^*s~^ 
corresponding to a wave number of around 1700cm~^. 
Hence, at room temperature, this oscillation is essentially 
in the ground state. This zero point motion gives rise to 
a mean square deviation of the position of the oxygen 
atoms of the order of 0.02A. The oxygen atoms in the car- 
bonyl group carry an effective negative charge and hence 
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their motion will affect the potential that charges such as 
potassium ions experience in the ion channel. The fluc- 
tuating potential leads to fluctuating energy differences 
between two potential wells which in turn are responsible 
for dephasing of superposition states between two wells. 
The potassium ions in their potential wells have a fre- 
quency of around ljk ~ 3 • 10^^s~^. Numerically one can 
now determine that the zero-point motion of the carbonyl 
groups will affect the energy difference between two wells 
by around fkoK/W for a typical displacement. This leads 
to a coupling of the oscillations of the carbonyl groups 
to that of the potassium ions. Now we treat this as a 
simple two-mode quantum system, one mode referring 
to the oscillation of the carbonyl group and the other 
to the energy difference between two potential wells. The 
Hamiltonian is given hj H = huicoci' a+huJK{Xco)b^b = 
hucoo-^ + 0.1/ictJ_K-(ffl + a^)h''h + hujKb^b. Here hujK de- 
notes the energy difference between two potential wells 
which is in turn a function of the position of the oxy- 
gen atom in the carbonyl group. Now initialising the 
system in the state |0)a(|0) -|- |l))i, corresponding to a 
superposition of an ion in the left or right potential 
well we find that the time in which the fidelity with 
the initial state has dropped to 1/e is around 3000a;^^ 
suggesting that the decoherence time is of the order of 
tdec ~ SOOOcj^^ ~ 3 • 10"^s. This is of the order of the 
tunneling time. It should be noted that there are a va- 
riety of unknowns that may affect the decoherence time. 
On the one hand we have ignored the oscillation of the 
carbon, nitrogen and hydrogen in the ion channel which 
may lead to a slightly higher decoherence rate. We have 
neglected oscillations along the channel axis as these have 
a much smaller effect on the potential and hence couple 
more weakly to the potassium ions. On the other hand 
the present estimate assumes a first order effect of the 
oscillation of the carbonyl group on the potential in the 
ion channel. If the system was making oscillations around 
equilibrium the potential change would be a second order 
effect and hence smaller. Hence the quoted figure provides 
a ball park figure but may vary somewhat. Furthermore, 
it should be noted that these oscillations do not only lead 
to a dephasing between two potential wells but also af- 
fects the height of the potential barrier. In fact, the line 
width A(jj due to the oscillations and the variations in the 
height of the potential barrier are approximately propor- 
tional. As the tunneling rate is strongly non-linear with 
the potential barrier height it can be expected that these 
oscillations can actually be beneficial. 
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